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Abstract: We have proposed an approach to the interference phase extraction in the 
homodyne laser interferometry. The method employs a series of computational steps 
to reconstruct the signals for quadrature detection from an interference signal from a 
non-polarising interferometer sampled by a simple photodetector. The complexity trade-off 
is the use of laser beam with frequency modulation capability. It is analytically derived 
and its validity and performance is experimentally verified. The method has proven to be a 
feasible alternative for the traditional homodyne detection since it performs with comparable 
accuracy, especially where the optical setup complexity is principal issue and the modulation 
of laser beam is not a heavy burden (e.g., in multi-axis sensor or laser diode based systems). 

Keywords: laser interferometry; optical metrology; homodyne detection; digital 
signal processing 



1. Introduction 

The laser interferometry represents the most precise technique for measurement of geometrical 
quantities today. Among currently known methods, it provides the best resolution and accuracy and 
ultimate (theoretically an infinite) dynamic range. The basic resolution of the laser interferometer 
corresponds to the wavelength of the laser source, but for many applications such an accuracy is 
insufficient. The resolution is further improved by various methods, commonly referred to as fringe 
subdivision or detection techniques [1,2]. 
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Generally speaking, the detection methods usually expand the basic concept of interferometric 
measurement by adding extra information to the interference signal that helps to extract the phase 
information with an increased precision [3,4]. A typical example is the homodyne detection that handles 
the polarisation of the laser beam in the optical hardware so that a compound wave is produced instead 
of the interference signal. The phase information is then extracted using additional electro-optic system 
that produces a pair of interference signals in quadrature (see details in the next section). The phase is 
then extracted from these signals with significantly increased precision. 

In this paper, we present a novel detection technique, related to a traditional homodyne detection, 
that relies on a digital computation instead of rather complex optical processing. It employs a frequency 
modulated laser source on the input of a non-polarising laser interferometer and a simple photodetector 
that observes the output. Then a combination of synchronous demodulation, linear transforms, mixing, 
phase shifting and trigonometric transforms is used to reconstruct the same pair of quadrature signals. 
The detection is further improved by a scale linearization techniques that mitigate the measurement error 
introduced when the method is applied in a real measurement system. 

The presented approach brings several advantages, especially reduced demands on the precise 
construction, reliability and cost of laser interferometry based measurement systems. Our aim was not 
to present a novel technique that would allow for better precision, accuracy or bandwidth — especially 
the bandwidth is limited by the frequency modulation bandwidth of the laser source (we deal with 
DC— 100 Hz in our experimentation). Instead, we proposed a method that achieves a performance 
comparable to that of the current methods, i.e., decent resolution (sub-nanometre). 

The rest of this paper is organised as follows. Section 2 summarises the principles of the homodyne 
detection and presents an analytic description of the proposed method. Section 3 describes the 
experiments we have carried out to validate the performance of our method in comparison with a 
homodyne detection, while Section 4 presents the results and holds the discussion. Finally, Section 5 
summarises and concludes the paper. 

2. Methods 

2.1. Homodyne Detection with Polarising Optics and Quadrature Detection 

In a typical Michelson setup with homodyne detection, the interferometer setup contains a polarising 
beam splitter to produce the measurement and reference beams with a different polarisation plane. As a 
result, there is a complex light beam compound from the waves with mutually perpendicular polarisation 
on the output of the interferometer rather than a simple interference. That compound beam is further 
processed by an optoelectronic detection system, usually referred to as a quadrature detection unit 
(schematically shown in Figure 1), that produces a pair of interference signals I x and I y that are in 
quadrature, i.e., mutually phase shifted by n/2. 
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Figure 1. Polarising laser interferometer with quadrature detection unit for homodyne 
detection: the interferometer employs the polarising divider, the reference corner cube 
reflector CC re / and measurement reflector CC meas . The output beam is compound from two 
beams with a different polarisation. The detection unit splits the beam into two detection 
branches on the non-polarising divider NP. In the first branch, the polarising beam splitter 
PBS 1 adjusts the polarisation to a single plane so that the interference occurs and is observed 
by photodetectors Dl, Dl' (mutually in an opposite phase). In the second branch, one of the 
wave components is retarded by the retarder plate RP. The photodetectors D2, D2' observe 
the interference that occurs due to rotation of the polarisation plane on the polarising beam 
splitter PBS2. The subtractor SI produces the signal I x as a difference between signals Dl 
and Dl', the subtractor S2 produces I y from D2, D2' in the same way. 



Polarising divider Laser 




The immediate value of the interference phase <p{t) is extracted by a quadrature detection, i.e., the 
interference phase is calculated using an inverse function 

#)=arctan(||) (1) 

from immediate values of the X- and Y-axis signals I x and I y . 

One of our objectives was to propose a method that would be simpler in terms of optical hardware 
complexity. This complexity is represented by, e.g., the polarising beam splitters, the corner reflectors, 
semi-transparent mirrors and the demands for optical system adjustment. This is why using exclusively 
non-polarising interferometer with planar reflectors and a simple photodetector can be considered 
interesting direction towards cost-effective measurement systems, despite the necessity of laser beam 
modulation. Several extensions of the basic scheme of the homodyne detection are currently used in 
modern measurement systems, e.g., the use of planar mirrors [5] or advanced retroreflectors (Mi-series 
interferometer, SIOS Messtechnik) instead of cube corner reflectors or the multichannel photodetector 
capable of producing the quadrature signals (Renishaw). 

The next section presents the analytical description of the interference phenomenon and the signals 
produced by the quadrature detection unit. Note that the laser beam is considered as a set of 
monochromatic plane waves in vacuum since an ideal source of coherent light is assumed. 
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2.2. Quadrature Signals 

The electromagnetic wave generated by a laser source can be described at a given point by a vector of 
electric field E that oscillates in the plane perpendicular to the direction of propagation. If the direction 
of oscillations is fixed we call the wave linearly polarised and we can describe it by the magnitude of 
electric field E. For definiteness, assume that the laser source is placed at the origin of the Cartesian 
coordinate system and the wave propagates along the x axis. 

The electric field at the source is given by the real part of the complex function 

E(t) = E 0 expli^t)] (2) 

where E 0 is amplitude and <fi(t) is a general phase. By definition, the phase is related to the instantaneous 
angular frequency u(t) by 

Hence the phase can be calculated from the known angular frequency as 

<j>{t) = J u(t)dt + 0o (4) 

The propagation of the wave from the source to the point x = I takes the time r = l/c, where c is a 
velocity of propagation. Therefore at the point x = I the electric field is 

E(T,t) = E 0 exv[i<f)(t-T)) (5) 

where 4>(t — r) is a phase shifted by time interval r. 

The intensity observed by a detector placed at the point x = I is proportional to the squared absolute 
value of electric field and is given by 

I(r,t) = q\E(r,t)\ 2 = qE(r, t)E*(r, t) (6) 

where the asterisk denotes the complex conjugation and q is a proportionality constant. 

Now assume there are two waves propagating along a single axis (practically representing two waves 
that travel along two arms of the interferometer). At the output of the interferometer the two waves are 
recombined, the interference occurs and the compound wave, usually denoted as an interference signal, 
is observed. Let l\ be the distance travelled from the source to the detector by the wave that propagates 
through the reference arm and l 2 the distance through the measurement arm. Then the electric field of 
the compound wave at the detector is 

Ei(T!, t 2 , t) = Ex(t x , t) + E 2 (t 2 , t) = E Q1 exp[i0(t - n )\ + E 02 exp[i0(t - r 2 )] (7) 

where T\ = h/c and r 2 = l 2 /c. 

The intensity at the point of observation can be calculated using Equation (6). Denoting I x = qE^ 
and I 2 = qE^ 2 and using the Euler's identity for complex numbers we obtain the observed intensity 



/(n, r 2 , t) = h + I 2 + 2JIJ 2 cos[0(t - n) - <f>(t - r 2 )] (8) 
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Changing the variables t — T\ — > t and denoting r 2 — T\ = r, we rewrite the observed intensity as 

I{r, t) = h + I 2 + 2^/hhcos[(f){t) - 4>{t - r)] (9) 

We see that the observed intensity depends only on amplitudes of two interfering waves and their 
phase difference. The latter is caused by the path difference I — l 2 — l\ between the paths travelled by 
two interfering waves. 

When the laser source works on single frequency /, the angular frequency of the wave is constant and 
can be calculated as 

to = 2nf = 2tt^ (10) 
A 

where A is wavelength. From Equation (4) we get the phase of the first wave 

<f>(t) = out + 0o (11) 
where 0 O is the initial phase. Similarly, the shifted phase of the second wave is 

(f>(t-T) = L0(t - T) + <f) 0 (12) 

The intensity I x observed by the first branch of the detection unit can be calculated by means of 
Equation (9). We get 

4(r, t) = h + I 2 + 2yflJ 2 cos(uT) (13) 

which is a well-known result that allows one to observe interference fringes when the time delay 
r = (I2 — h)/c is not constant. This can be easily achieved by changing the length l 2 of the 
measurement arm. 

The intensity I y observed by the second branch of the detection unit (the Y-signal) that involves the 
retarder plate is shifted from I x by the phase factor n/2 and is given by 

I y (r,t) =I 1 + I 2 + 2^LJ 2 sm{uT) (14) 

On the basis of analytic model of the quadrature signals, the next subsection presents the approach to 
computational reconstruction of quadrature signals from a signal observed by simple detector. 

2.3. Reconstruction of Quadrature Signals 

As mentioned above, we employ a frequency modulated laser beam instead of a monochromatic 
one. Similar approach has been already employed, e.g., in phase-shifting and phase-modulation based 
interferometry [6-8], frequency modulated continuous wave interferometry [9] or for absolute distance 
measurement [4,10], nonetheless the follow-up evaluation is principally different from our approach. 

Particularly, we apply the sinusoidal frequency modulation. First, we denote the mean value of 
angular frequency of the wave cu s and the amplitude of the sinusoidal modulation cu m . Then the angular 
frequency of the first wave is 

u(t) = u s + u m sin(Q£) (15) 
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where Vt is modulation angular frequency. The phase of the first wave can be obtained by integration of 
the angular frequency. We get 

0(t) = J u(t)dt + 0 O = u s t - ^ cos(fit) + 0o (16) 

Similarly, the phase of the second wave shifted by the time delay r is 

0(t - r) = u a (t - t) - ^ cos(fi(t - r)) + 0 O (17) 

The intensity I x observed by the first branch of the detection unit can be calculated using Equation (9). 
For clarity, we first calculate the phase difference 0(i) — 0(i — r). Using the standard relations for 
trigonometric functions we obtain 



0(t) - (p(t - r) = u s r + sin 



■ At - T 2 



.VLt\ 

sin — (18) 



2 J 

This expression can be further simplified on the basis of the following considerations. Due to the 
typical dimensions of commonly used interferometers and a large value of the speed of light, the time 
delay r is negligible compared to the time t. Moreover, the term Vtr <C 1. Hence we can make the 
approximations 



sm 



<Ht- T - 



'Q T \ fir 

sin(Of). sin(— J=— (19) 



The phase difference is then 

0(t) — 0(t — r) = uj s t + u m T sin(fit) (20) 

Inserting this expression into Equation (9) we finally get the intensity I x observed by the first branch 
of the detection unit 

I x (r,t) = h + I 2 + 2yJTj~2 cos [u s t + u m r sin(fii)] (21) 

The intensity I y observed by the second branch of the detection unit differs from I x only by global 
phase factor n/2, therefore we have 



I y (r,t) = h + h + 2\jlxh sin[u; s r + u m r sin(fit)] (22) 

We see that the observed intensity is variable even if the time delay r is constant. This is obviously 
caused by the modulation of the angular frequency of the wave emitted by the source. 

The intensities I x and I y observed by particular branches of detection unit given by Equations (21) 
and (22) were derived for the general case of sinusoidal modulation of the angular frequency. In these 
calculations the depth of modulation tu rn and time delay r were not subjected to any restrictions. 
However, we show that under certain assumptions for these two variables, the Equations (21) and (22) 
can be substantially simplified. 

Since the absolute value I\ + I 2 and the factor 2^I\l2 are common for both I x and I y , we omit them 
for better readability of the following text. Mathematically, this corresponds to the normalisation of 
functions I x and I y by means of relations 
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Then we start with normalised intensities 

I x (t, t) = cos[u; s t + uj m T sin (fit)] 

I y (r,t) = sin[u; s r + u m r sin(fit)] (24) 

Using standard formulae for trigonometric functions we can rewrite them to 

I x (t, t) = cos(u; s r) cos[u m r sin(fit)] — sin(co><,r) sin[a; m T sin(fit)] 

I y (T,t) = sm(u s r) cos[u m r sin(fit)] + cos(w s r) sin[a; m T sin (fit)] (25) 

Since the maximum absolute value of the term sin (fit) equals to unity, assuming that u m T C 1 we 
can make the following approximations 

cos[co> m r sin(fit)] = 1, 

sin[u; m T sin(fit)] = u; m r sin(fit) (26) 
and for observed intensities we get formulae 

I x (r,t) = cos(u s t) — u m r sin(fit) sin(o; s r) 

I y (r,t) = sin(a; s r) + u m r sin(fit) cos(a; s r) (27) 

Both observed intensities I x and I y are now expressed as a sum of two terms. We explain the meaning 
of particular terms for intensity I x . The discussion for I y can be done in a similar way. 

The first term of I x expressed as cos(o; s r) depends only on the mean value of angular frequency 
uj s and time delay r. Hence it is identical to the corresponding term calculated for interference of 
two monochromatic waves, see Equation (13). The second term of I x is the product of amplitude 
co m Tsm(co s r) and oscillating function sin(fit), which corresponds to the modulation of angular 
frequency. Since we have assumed that co m r 1, the amplitude of the second term of I x is much 
smaller than the amplitude of the first term. 

For further considerations we make the following assumption. When the mirror in measurement arm 
is moving we demand that the function sin(fit) oscillates much faster than cos(a; s r). Mathematically, 
this corresponds to the relation 

d(fit) d(u s T) 

-V > -ir (28) 

For illustration, considering the motion of the mirror in measurement arm with a constant speed v, 
the last condition is expressed as fi 3> co s v/c. 

The mean values (I x ) and (I y ) of observed intensities can be calculated using the relations 

(I x (r,t)) = — I x ( T ,t)dt = - cos(u; s t) - w TO rsin(fit) sin(uvr) dt 
1 Jo 1 Jo 

1 r T ~ 1 r T 

(Iy{T,t)) = — I y (r,t)dt = — [sin(w s r) + w m rsin(fit) cos(w s r)]dt (29) 

where T = 27r/fi is a period of frequency modulation. The stated assumption Equation (28) means that 
the value of functions sin(w s r) and cos(uj s t) can be considered constant in time interval from 0 to T. 
Moreover, since the integral of sine and cosine function over the period is zero, we get the mean values 

(4(r,t)) = cos(w s r), 

(I y (r,t)} = sin(w s r) (30) 
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We see that under the stated assumptions, the mean values of observed intensities when the 
frequency is modulated are identical to the intensities observed when the frequency is constant, see 
Equations (13) and (14). Note that the suitable depth of modulation is crucial for the reconstruction of 
the quadrature signals. 

Now we would like to show the procedure that allows us to calculate the mean value of intensity (l y ) 
from the observed intensity I x . We start with intensity I x at the time t given by 

Ix(T,t) = cos(u; s t) — u m r sm(ftt) sin(a; s r) (31) 

and its value at neighbouring time t — At, which is 

I x (t, t — At) = cos(u s t) — u m T sin(0(t — At)) sin(w s r) (32) 

Calculation of their difference AI X gives us 

A4(r,t) = I x (r,t) - I x (r,t - At) = 

= u m T sm(u s T)[sin(Q(t — At)) — sin(Qt)] = 

= — u m T sin(o; s r)[(l — cos(f2At)) sin(f2t) + sin(f2At) cos(f2t)] (33) 

Next, multiplying the difference AI X by cos(f2t) (note that it is a phase-shifted copy of modulation 
signal, originally sin(fit)) and calculating the mean value we get the derived signal 

I d (r,t) = (Ai x (T,t)coa(Qt)) = 

= -^^sm(uj s T) [ [(1 - cos(fiAt)) sin(fii) cos(fit) + sin(fiAt) cos 2 (Qt)]dt 
T J o 

= —-u m T sin(QAt) sm(u s r) (34) 
Finally, from comparison with Equation (30) the derived signal I d is 

I d (T,t) = --u m Tsm(nAt)(I v ) (35) 

We see that the derived signal Id is proportional to the mean value of intensity (I y ). Hence, we are 
able to reconstruct the mean intensity (I y ) from intensity I x . 

The last part of the novel detection method, presented in the final part of this section, addresses the 
issues with the practical implementation. 

2.4. Scale Linearity Compensation 

We have identified several sources of inaccuracies related to a practical implementation of the 
proposed method: 

(i) the measured intensities I x and I y are generally burdened with a measurement error, e.g., noise or 
the quantisation error (i.e., the error caused by the discretization of analog measure), 

(ii) the currently available laser sources are liable to a residual amplitude modulation (RAM) [11], i.e., 
the modulation of their optical frequency has an unwanted influence on the intensity and spatial 
characteristics, 



Sensors 2012, 12 



14103 



(iii) the computation over discrete representation of the analog signals causes additional inaccuracy, 
e.g., numeric calculation of derivation, round-off error. 

Note that the list of influences is far from being complete, nonetheless these were the most significant 
ones that we met and found important to be addressed. 

Currently, the problems with scale linearity (i) are well explored [12,13] and there are several methods 
to deal with them [14-16]. 

The RAM effect (ii) is superimposed onto the effect of the laser frequency modulation. As a 
consequence, there occurs an additional phase shift between I x and I y besides the desired ir/2 when 
our detection method is used. This results in an unwanted offset in the detected interference phase that 
depends on the difference of path-lengths and the modulation depth. The effect is additionally amplified 
due to the discrete computation (iii), especially due to the calculation of the numeric difference over a 
discrete set of sampled points. 

The phase shift is typically corrected by the scale linearization techniques, nonetheless our 
preliminary experimentation indicated that a compensation of the phase shift between the quadrature 
signals before the elliptic-fitting based linearization technique leads to better overall performance of 
our method. On the basis of this experience, we have employed both techniques: the phase shift 
removal is described in the following section and for the scale linearization technique we used an existing 
method [14] — the latter method is not described in this article. 

2.5. Removal of Unwanted Phase Shift between Quadrature Signals 

The principle of our phase shift removal method is the following: after calculation of I y from the 
observed I x we detect their mutual phase shift 5 a (caused mainly by RAM), then we transform I y in an 
appropriate manner to correct this phase shift and after that we use the linearization method. 

The following description presents the method for a general case, i.e., we assume unequal amplitudes 
of the quadrature signals, even though we have operated with unit amplitudes of I x and I y before. 

We start again with intensities I x and I y in the form 

I x = K x cos 0, 

I y = K y sin(0 + 5 a ) = K y cos + 5 a - 0 (36) 

where we assume that their mean value is zero and K x , K y are real amplitudes of I x , I y . Note that the 
amplitudes are known from the observed time course of I x and I y in practical applications. Next, we 
introduce a quantity K xy that corresponds to the amplitude of the sum I x + I y , formally written as 

K xy = m&x(I x + Iy). (37) 

The value K xy is also known from observed data. 

We illustrate the functions I x and I y given by Equation (36) in a phasor diagram using two phasors 
K x and K y of lengths K x and K y , respectively (see Figure 2(a)). These phasors together form a constant 
angle 7r/2 — 5 a resulting from Equation (36), but they both rotate around the origin as the phase <p 
changes. In such a diagram the values I x and I y are simply given by orthogonal projection of the phasors 
K x and K y to the x-axis, respectively. 
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Figure 2. Detection of the unwanted phase shift 8 a : the known instantaneous values I x , I y 
are projected to phasors K x , K y that form the angles 0 (equal to the immediate interference 
phase) and 0' = 0 — it/2 + 5 a with axis x (a); a combination of K x and K y forms phasor 
(b). Since there exists such 0 such that K xy is parallel to axis x, the real angle 
a = 7r/2 + 5 a can be calculated using law of cosines. 




The geometrical meaning of quantity K xy is the following: since it was defined as the amplitude of 
sum I x + I y , it represents a length of phasor K xy = K x + K y in the phasor diagram, see Figure 2(b). 
To derive the expression between the lengths K x , K y , K xy and phase difference 5 a , we use the law of 
cosines. Since the angle a shown in Figure 2(b) equals to 5 a + tc/2, we get 



K xy = K 2 X + K y - 2K x K y cos [5 a + - J (38) 

It should be emphasised that the function K xy (S a ) is monotonic in the range — n/2 < 5 a < n/2 for 
arbitrary values of K x and K y , which can be easily verified by calculating the derivative of K xy with 
respect to 5 a . This allows us to solve Equation (38) with respect to 8 a . We obtain 

(K 2 -K 2 -K 2 \ 
5a = arCSin { XV 2K x K y ~ ) (39) 

Since the value K xy and amplitudes K x , K y can be easily determined from measured data, we are able 
to calculate the phase shift 5 a and eliminate the unwanted phase shift between the quadrature signals. 

For this purpose we now need to transform intensities I x , I y given by Equation (36) to another pair 
4> r y g iven b y 

I' x = K x cos 0, 

/; = 7^sin0. (40) 

Here I' x is chosen to be identical to I x , for I' y we need to derive the transformation relation expressed 
in terms of I x , I y and 5 a . Therefore we rewrite Equation (36) as 

I x = K x cos 0, 

K 

I y = K y (sin 0 cos 5 a + cos 0 sin 5 a ) — I' y cos S a + I x —^- sin 5 a (41) 

■"■X 
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Solving the second equation with respect to /' 



we obtain the transformation relations 



I XI 

I V K X - I x K y sin 5, 




a 



(42) 



The transformed intensities I' x and I' are now deprived from the additional phase shift 5 a caused by 
amplitude instability of the laser source. As a next step, the signals are scaled to unit amplitude, the scale 
linearization is applied and then the phase is extracted. 

Note that the algorithm is principally based on a single data-point determined from the measured data 
(maximum value of I x , I y ,...\ so that some outlying data points might cause a substantial error in the 
calculations. Since the K x and K y (and 5 a as well) can be considered constant, the method can rely on 
some kind of statistical information (and actually has within our experiments). 

3. Experimental Verification 

We have carried out several experiments to verify the performance of the proposed technique. 
We have assembled an optical setup that allowed us to make the displacement measurements at the 
nanometre scale referenced to an interferometer with the homodyne detection technique, i.e., quadrature 
detection unit. 

The reference was physically provided by a calibrated quadrature detection unit in combination 
with the support logic. The reference delivers the reference phase information with measurement 
uncertainty of AT/1, 024 (where iV denotes single interference fringe), i.e., 0.35° [17] (including the 
linearity compensation). 

3.1. Experimental Setup 

The experimental setup, schematically shown in Figure 3, comprised a Michelson interferometer with 
polarising optics and two detection branches: the reference branch with a quadrature detection unit and 
the testing branch that performed the proposed technique for interference phase detection. 

A solid-state laser (frequency doubled Nd:YAG at 532 nm with linearly polarised beam; 
PROMETHEUS by Innolight, Germany) was used. The output beam was collimated to particular 
diameter (3 mm) that optimally fits the detection unit and the photodetector. A half- wave plate (A/2) 
was used to adjust the polarisation plane. The measuring reflector was mounted on a piezo actuator and 
on a motorised stage. This combination allowed for the measurement on an operational range of 600 mm. 
Note that the polarising variant of the beam splitter is used, as it is a necessary condition for the reference 
interferometer (the detection apparatus, not to be confused with the reference arm of the interferometer). 
The proposed method does not require preserving different polarisation in individual arms. 

A non-polarising splitter (NP) is used to form two detection branches: the testing branch that involves 
our novel detection method with an ordinary photodetector (PD) and the reference one that uses the 
traditional detection unit (denoted QDU). 
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Figure 3. Experimental setup scheme: the interferometer employs the polarising divider, 
the reference corner cube CC re / and measurement corner cube CC meas . The laser input is 
attenuated by a neutral density absorption filters ND, collimated by an optical expander EXP 
and its polarisation plane is rotated by the half-wave plate A/2. The interferometer output 
is split into two detection branches by the non-polarising divider NP. In the first branch, the 
interference occurs on the polariser POL and the intensity is observed by a photodetector 
PD. In the second branch, the signal from the interferometer is processed by the reference 
quadrature detection unit QDU. The acquisition and control hardware drives the modulation 
of the laser source optical frequency and evaluates the signals from both the photodetector 
and the quadrature detection unit. 
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The bandwidth of both detection methods is limited by the combination of several factors: response of 
the photodetector, signal amplifiers, analog-to-digital conversion sample rate and processing capability. 
In our particular case, the available bandwidth was 50 kHz for the QDU and 100 Hz for the novel method 
{i.e., 100 phase readings by QDU and 200 phase readings by tested method per second). The optical 
line-up in the testing branch is extended by a polariser plate that alters the orthogonal polarisations into 
single plane so that testing branch observes the interference signal instead of compound wave (that is 
processed by the reference branch). The bandwidth of our novel method is relatively low, nonetheless 
the sense was to demonstrate the principle and compare the linearity with a reference method. 

The signal generation (for laser source modulation) and acquisition was done by a dedicated 
DAQ device (U2531 by Agilent). The generation part produced synchronised sine and cosine wave — the 
sine controlled the modulation, the cosine was synchronously sampled (acquired) together with the 
output of the QDU and the PD output and used for the quadrature signal reconstruction. The reason 
to generate the modulation signal and its phase shifted copy and then sample the copy again was to 
achieve precise synchronisation between the generated outputs and signal inputs. 

3.2. Design of Experiments 

The measurement arm hardware was successively placed into different positions (measurement 
points) 

AL e {0, 1, 2.5, 5, 7.5, 10, 15, 25, 37.5, 50, 62.5, 75, 100, 150, 200, 250, 300, 400, 500, 600} mm 
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to cover the operating range AL e [0, 600] mm. At each point, the measurement arm was continuously 
displaced (with approximately constant velocity using the piezo actuator) over the trajectory of /im 
(more than two complete fringes) and the detection outputs were acquired. The selection of points 
was suggested by preliminary experiments that indicated that the correct distance determination 
is more problematic at the beginning of the operating range. The zero position (AL = 0) 
actually denotes a closest achievable position of the measurement arm, given by the mechanical 
construction of interferometer components and corresponds actually to the path length difference of 
approximately 2 mm. The laser source has been modulated by a free running sine wave (at 1kHz). 
The depth of modulation was empirically estimated and manually configured for each point — the optical 
frequency was tuned approximately by ±20. ..50 MHz. 

The outputs from the quadrature detection unit I x , I y , the photodetector signal I p and modulation 
signal M in were synchronously acquired (at sample rate of 100 kHz). One phase reading of the 
novel method took a time period that corresponds to five complete modulation periods (5 ms). Three 
measurement cycles were recorded at each measurement point. The evaluation was performed afterwards 
over the entire cycles at once. 

We have inspected the difference between the interference phase information extracted by the 
reference and by the testing method, referred to as the phase determination error (PDE). In the reference 
branch, the X- and Y-axis signals from the quadrature detection unit were recorded and the scale 
linearization was applied. The photodetector output I p was used for reconstruction (standing for the I x 
when related to the previous analytic description) of the quadrature signals using the presented method. 

4. Results and Discussion 

4.1. Results 

A total of 60 measurement cycles in 20 positions was completed, where each cycle contains 890 phase 
readings. A sample evaluation of a single cycle is shown in Figure 4. 

There is an apparently good coincidence between the quadrature signals from reference and testing 
detection branch (compare Figure 4(a) to Figure 4(b)) (the slight fluctuations are caused by changes 
in refractive index of air, but since the optical paths are shared, the influence is subtracted within the 
comparison). The phase extraction error plot (see Figure 4(c)) indicates that the PDE keeps within ±1° 
range (for the particular sample). 

The linearity of our method is summarised in Figure 5. There is an observable periodic error, and 
the mean PDE fluctuates within ±0.15° . The Figure 6 shows the overall dependence of the mean PDE 
on the path length difference when the results from all three cycles measured at the same distance point 
were combined so that the repeatability aspect is incorporated. The overall mean PDE stays below ±0.5° 
and the standard deviation stays below ±0.75°. The jitter of the mean error over the operational distance 
can be attributed to the optical adjustment. The combination caused a slightly greater scatter of the phase 
determination error. 

Several related issues were omitted within the experimentation, e.g., the refraction index of air 
influence [18], the uncertainty of the laser source stability or a precise temperature stabilisation. In 
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our case, when two detector paths were compared, these influences were suppressed by the differential 
measurement, since the optical paths and signal acquisition within the detection chain were common for 
both methods. The results show that the novel method has a comparable accuracy to the reference 
method. Exploring the limits of our proposed method is still a subject to a future research and a 
comparison with a more accurate method is required to state the absolute precision and uncertainty. 

Figure 4. Evaluation of single experimental cycle: the X- and Y-axis signals from quadrature 
detection unit, i.e., the (I x (r,t)) and (Iy(r,t)) (refer to Equation (30)), with scale linearity 
compensation applied (a); X- and Y-axis signals reconstructed by proposed method, i.e., 
(I x (t, t)) and Id (Equation (35)), reconstructed from the signal observed by the photodetector 
with the phase shift removal and scale linearization applied (b) and the phase determination 
error, PDE (c). 



a) 



b) 



c) 




Figure 5. Experimental results: the red points show the periodic error, i.e., the dependence 
of the mean phase determination error on immediate interference phase; the blue error 
bars indicate the interval of radius of standard deviation around the mean PDE at each 
corresponding point i (n* ~ 1, 485) 
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Figure 6. Experimental results: the red points show the dependence of the mean phase 
determination error on the path length difference; the blue error bars indicate the interval 
of radius of standard deviation around the mean PDE at each corresponding point i (rii ~ 
2,670) 
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4.2. Discussion 



The results indicate that the proposed method is a suitable alternative to the traditional homodyne 
detection, especially in situation when the accuracy is not the most critical aim but the construction 
simplicity and scalability is. There are several issues that need to be considered to judge the novel 
approach: its dynamic range, the trade-off between optical complexity and the requirement for laser 
modulation, the response bandwidth limitations and the uncertainty introduced by the computations. 

As the homodyne detection has theoretically infinite dynamic range, it allows for the measurement 
with the same precision even on long distances (several kilometres). The practical limitation is given 
only by the construction issues, laser stability and influences of the environment. The novel method is 
principally limited by the available depth of laser frequency modulation [19], but using an approach for 
an adaptive control of the modulation depth the method would be probably applicable on a comparable 
operational range. 

There are also higher demands posed on the laser source — there is the need for the laser beam 
frequency modulation. The modulation bandwidth limits the response of the measurement systems since 
the duration of at least one complete modulation period (five in case of our experiments) is processed to 
produce a single phase reading. Consequently, velocity of measurement is significantly slower than the 
optical processing. The demand for laser with a fast modulation capability can be fulfilled, to a certain 
extent, by cost-effective semiconductor lasers [20,21]. 

On the other hand, the method requires less complex optical setup to achieve comparable 
results — there is no need for polarising beam splitters and semi-transparent mirrors within the optical 
lineup. As a consequence, the method also scales well, e.g., when several measurement axes would 
be required within a measurement system. In such a case, the added complexity per each added axis 
is minimal in comparison to the traditional means — it has been already demonstrated in a practical 
application for surface diagnostics [22] and multi-axis geometric measurements [23]. An extra axis 
requires only a small amount of additional optical components and poses no substantial additional 
requirement on the laser source (assuming there is enough power available, but it is usually not an 
issue, e.g., for green Nd:YAG lasers). 

In our experiments, the modulation bandwidth was limited by available high-voltage amplifiers to 
200 phase readings per seconds (corresponds to 1kHz modulation). With the particular laser source, 
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a 100 x higher response would be achievable; further enhancement would be possible using and 
acousto-optic modulator. 

In the basic principle, similar computational methods are also sensitive to the noise issues. In case 
of our novel method, the issue is well mitigated by a massive averaging and suitable signal filtering 
without increased demands on system performance as the operations are essential parts of the phase 
extraction method. 

There were also several assumptions we stated (in sake of simplicity) in the physical description of 
the novel approach. We assume the time delay r — the difference in wave propagation along unequal 
optical paths in the arms of the interferometer — is negligible (see the paragraph above Equation (19)). 
In the borderline case, the path length difference 1 = 1, 200 mm, so the time delay r « 4 x 10~ 9 and 
the term fir k4x 10 -6 . We also assume that the product of r and the modulation depth co m r <C 1 
(Equation (26)). In the case when I = 1, 200 mm and co m = 20 MHz, the product co m T = 0.08, which 
corresponds to the limit of validity of our assumption co m r <C 1. Finally, we assume that the phase 
change (due to displacement) is slower than the modulation of the laser, i.e., fi ^> co s v/c (Equation (28)). 
The displacement velocity v was v ~ 2.5 x 10~ 7 m/s, the Q = 1kHz, thus for the mean frequency 
co s « 563.25 THz the ratio is 1 x 10 3 > 7.5 x 10^ 6 . 

5. Conclusions 

We have proposed and verified a novel detection technique that represents a novel approach to the 
homodyne detection in laser interferometry. 

A detailed analytic description of the principles has been presented and the validity and performance 
were evaluated by experimental means. 

The experimental results indicate that the novel approach performs comparably to the traditional 
detection technique. Since it requires simpler optical hardware setup — e.g., exclusively a non-polarising 
interferometer, single photodetector, but frequency modulated laser beam and computational power — it 
can be considered a suitable alternative to the techniques that rely on additional optical processing. 
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